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A description of the thermodynamic properties in the critical region <>( a physical system is ob- 
tained from a scaled expression lor the free-energy /(f). T). In jieneral, a mmsymmetric coexistence 
curve is predicted, with the symmetric case (e.ji., magnets) included as a special example. Kor fluids, 
deviations from symmetry j^ive rise to an expression for the average densit) below the critical point 
nonlinear in the temperature near T, (in contrast to the usual "law <>! rectilinear diameter"): these 
asymmetries also contribute to the discontinuity in the specific heal along the critical isochore. To 
lowest order, the formulation reduces to Widom's homogeneous scaling; the classical equations of 
state of the van der Waals type are incorporated as a special case. 
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1. Introduction 

The Widorfi-Kadanoff-Domb scaling hypothesis 
[1-3] ' provides a qualitative characterization for the 
thermodynamic features of a second-order phase 
transition [4]. Scaling incorporates the usual power 
law forms for the various physical anomalies and more 
generally serves as a lowest order asymptotic descrip- 
tion close to the critical point [5]. Such limiting be- 
havior has been found in various mathematical models 
and also appears to be valid in real systems. 

Despite the relative success of the scaling; approach, 
there remains a definite need for a more global picture 
of the critical region. This is evident during the analysis 
of bulk thermodynamic data where there exists first 
the problem of determining those physical quantities 
which satisfy the basic conditions required by simple 
scaling. For some systems, this selection is easily based 
upon some known/intrinsic physical property (sym- 
metry); whereas for others, the actual choice is less 
certain [5]. The determination is further complicated 
by the secondary problem, that of estimating the actual 
range over which scaling remains valid [6]. 

These difficulties are particularly evident in the 
application of the scaling description to real fluids near 
their critical point [4]. Here there is some experimental 
indication of a partial anti-symmetry about the critical 
isochore in the chemical potential p — number density 
p representation not present with either the pressure 
p or the specific volume v coordinates. It is also known, 
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however, that the phase boundary lacks this same 
symmetry in any set of variables. The coexistence 
curve is rather asymmetric about the critical isochore 
with such deviations being described by the empirical 
"law of rectilinear diameter." Despite the fact that 
this nonsymmetric behavior cannot be accommodated 
within the usual lowest order sealing description, the 
fluid data appear to scale in the preferred p,-p coordi- 
nates over a rather large region about the critical point 
(in a density range of ±30 percent and a temperature 
width of several percent from their reduced critical 
values [4]). 

An understanding and resolution of such behavior 
is hampered by an almost complete ignorance of what 
is beyond the lowest order scaling description. Since 
any results will be biased by the neglect of unknown 
terms, one cannot hope to extract any more than a 
simple qualitative picture from the experimental data 
without some knowledge as to the structure of correc- 
tions to the asymptotic scaling form. The present 
article represents an attempt to gain some insight into 
the structure of such terms. The idea was to construct 
a "model" calculation incorporating modification of the 
usual scaling formulation. Based upon an extension of 
classical thermodynamic ideas, the scheme describes 
systems lacking any known or simple symmetry and 
extends a scaling interpretation over an enlarged region 
about the transition. The approach allows for the first 
time discussion of the various contributions arising 
from asymmetry and higher ordered corrections to the 
asymptotic scaling behavior and their influence on 
observable physical properties. (After completing the 



103 



work reported here, we learned that Domb and Gaunt 
have looked at the structure of correction terms to the 
completely symmetric Ising model equation of state 
[71.) 

It should be pointed out, however, that the formula- 
tion given here is just one possible way of going beyond 
simple asymptotic scaling. For this reason, it is most 
important to emphasize the underlying formal structure 
rather than any of the individual numerical results. 
While the discussion is applied to fluid systems, the 
approach is easily adopted to other physical situations. 
The magnetic problem, with a completely symmetric 
phase boundary, has been discussed elsewhere [8]. 

2. Formulation 

Our approach is based upon two relatively simple 
observations. Just as in Widom's original hypothesis [1], 
classical thermodynamics suggests the framework of 
our analysis. First, recall that from a free-energy func- 
tion analytic in two physical variables, density p and 
temperature T, one will obtain an equation of state 
exhibiting in general a nonsymmetric phase boundary, 
deviations from a completely symmetric coexistence 
curve giving rise to a finite slope of the rectilinear 
diameter. With such behavior, however, one also ac- 
quires the classical (van der Waals-like) and quite in- 
correct description of the critical region [9J. Thus, 
instead of using such physical variables, we shall 
follow Widom fl] and assume the thermodynamic 
properties near the critical point are appropriately 
described in terms of one intensive physical quantity 
and one suitably reduced variable x. Moreover, the 
description is to be consistent with the homogeneity 
hypothesis, i.e., all functions are to be regular and of a 
general homogeneous form [5]. The existence of such 
nonclassical forms has been demonstrated in various 
model calculations where the scaling behavior appears 
as the first term in an expansion about the critical 
point. As an example, Cooper and Green [10] found 
that the appropriate free-energy of the Bose fluid is a 
function of the reduced temperature t= (T—T c )jT c 
and a reduced variable x=t | Ap | _1/ ^, Ap= {p — p c )lp c 
with the scaling form as the dominant term in an expan- 
sion containing higher ordered homogeneous terms. 
(Gaunt and Domb [7] have suggested a similar structure 
which applies to the Ising model.) As our second point, 
we shall assert that in general the free-energy density 
in the critical region is given by a convergent expansion 
of the form 



F(p, T) -F (p, 7 1 ) = I Ap | <* ' £ (Ap)-y,(x) (1) 



where x = t | Ap | _1/ ^. Here both 1//3 and 8 are to be 
arbitrary constants (not necessarily integer) and 
Fo(p, T) denotes any smooth analytic part of the free- 
energy near the transition p c , T c . It is implicitly as- 
sumed that the function F(p, T) is single value every- 
where and satisfies the well-known convexity conditions 
for thermodynamic stability [5]. The set of functions 



ifs(x)} are to be "well behaved" in the variable 
x=t | Ap | ~ ll(3 over the real interval of interest (min 
x { y) ^x^ + oo. (We require only that F(p, T) satisfy 
the physical conditions of continuity and differentia- 
bility in the single-phase region with each f s {x) 
analytic for x less than some finite constant R (see 
below).) Equivalently, one could take similar expres- 
sions for the chemical potential flip, T) and the pres- 
sure p(p, T) about the critical point p — p r , £ = 0: 

p(p,T)=p r + /x U)+ lApl 8 - 1 J (Ap) f /i.(*), (2) 



p(p, T) = Pc + po(t) + I Ap | s- 1 J (Ap)*p,-U). (3) 

.S'=l 

Both thermodynamic forces are to be well-behaved 
everywhere in the single-phase region and also on the 
phase boundary (except perhaps at the critical point). 
These expressions are also consistent with the experi- 
mental observation that both p and p have continuous 
first derivatives across the critical isochore. The domi- 
nant 5= 1 term of these expansions represent the usual 
asymptotic scaling expressions. The functions p s (x) 
and p s (x) are not independent but are themselves 
related to the set {f s (x)} through the thermodynamic 
relations 



p{p, T)=p 



, dF(p, T) 
dp 



p(p,T)= 1+p-j F(p. T). 



(4) 
(5) 



The form of Fo(p, T) must be consistent with the 
expressions for p and p. The first few important rela- 
tions of this hierarchy are that 



6+1 



x 3 
/3 dx _ 



fi(x)=p l (x) =p- l pi(x) 



p c l po{x)-p 2 (x) 



(3 dx _ 
2.1. Two-Phase Region 



Mx). 



(6) 



(7) 



Below the critical temperature, t < 0, we seek simul- 
taneous (non-trivial) solutions to the conditions for 
thermodynamic stability of coexisting phases, pi 
and p (; , 

p(p/,, T) =p(pGi T) andp(p/,, T) = p(p (; , T). (8) 
Denoting the two branches by p\— p c = Ap + and 



pa~ p r = Ap-, with Ap - = =t ( | r | 



3 , the shape of 



the phase boundary is found by expanding eqs (2) and 
(3) in a Taylor series about x ( f respectively, algebraic 
consistency requiring that 



x ± = x$ + xf\t\P-\-xi\t\ 2 P-\- . 



(9) 
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Matching the various powers in \t\& we find that 

(i) x$=Xo=— Xku with both p x {x) and P\{x) 

vanishing at x = — x () ; 
(ii) x( = — xt = Xi 

where x\ = I— Xq\ 1_/j 



PaU)— PdJLs(x) 



fi(8+l)-x—) 



(iii) in general #£ = (— l) r # r =x r > 

Thus we identify ^i(jr) with Griffith's [5] h{x) function 
and the condition that h(x = — Xo) vanish along the 
phase boundary. 

With the boundary condition fi\(—Xo)=0 we can 
also solve the differential equation (6) for f\(x) . The 
homogeneous solution is a constant times |*|^ (8+1) ; the 
general solution may be written as 



/,(*)=* w+l) 



/i(—*o) 



/8 



J-* V (6+ 



(y)_ 

i ) 



(10) 



The difference between the two coexisting densities 
is found to be given by 



pL-pG = 



2 



■x () 



+ 



PQB+1) tx x 

2 \#o 



-*2 
X() 



\t\»+. . 



(ID 



Deviations from a completely symmetric coexistence 
curve are described by the mean density of the two 
phases 



1 



P = 2 (Pl + Pg)=Pc + 



-/&i 



■x {) 



0+1 



|t|*+0(|*|*). 



(12) 



This is equivalent to the usual "law of rectilinear 
diameter" in the classical case (/3 =1 /2), but has the 
curious feature that for real fluids {/3 < V2) the average 
density p is curved and comes into the critical point 
with zero slope. See figure 1. When Xi = 0, the phase 
boundary is symmetric to this order. 

It is also possible to describe the form of the vapor 
pressure curve using the equation for p(p, T) 



Pcoex.(p, r)=Pc + p-(t) + |*| W+1) 



Pz(-Xq) 



to 



m+D 



+ 0(.^) 



(13) 



with a similar expression for the chemical potential in 
the two-phase region (replace p s by fi s ). Here p°(t) is 
a smooth analytic function in the temperature. 

While the specific heat C v in the two-phase region 
can be obtained from the free-energy F(p, T) , it is 




FlGl RE 1. Locus of mean densities for coexisting Phases. 

Curve I (dashed line) usual rectilinear diameter; 

Curve 2 (vertical line) line of mean densit) of completel) symmetric phase boundary 

like thai supposed by lowest order scaling; 
Curve 3 (solid line) locus of mean densities suggested by eq (12). 

perhaps more interesting to look at the individual con- 
tributions from the chemical potential /a and pressure 
p as exhibited in the thermodynamic identity (Yang 
and Yang 1 11|) 



iM-"(f£ 



(14) 



As the critical point is approached from below T c 
along the critical isochore p — p c , the anomalous 
temperature variation is given by 



PcC v 



A-- |d# 8+i >-2+o(f0< 8+3 >- 2 ; 



(15) 



where 

i4- = 0(8+l) {/3(S-f-l)~ 1} 



P'i(— xo)— PcP-ij— Xo) 



Xo 



]8(S+l) 



Here the expressions for the pressure (13) and chemical 
potential on the phase boundary have been used and 
it has been assumed that the difference p 2 U)— p*jM*) 
does not vanish atx = — x . The asymptotic behavior of 
the specific heat is thus characterized by a critical 
exponent a = 2-/3(8+ 1) and a coefficient determined 
entirely by the form of/i(jc) at x = — x {) (see eq (7)). In 
the simple scaling formulation where fx(p, T) is 
regular, the divergence of the specific heat along the 
critical isochore (T<T r ) comes entirely from the 
(d 2 pldT 2 ) Pc contribution. This formulation permits con- 
tributions from either the pressure and chemical 
potential by treating both thermodynamic forces on 
an equal basis. 
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2.2. Single-Phase Region 

All of the thermodynamic properties of the homoge- 
neous fluid in the one-phase region can be obtained 
directly from the free-energy F{p, T) . We assume that 
no higher order discontinuities exist anywhere in this 
region including the critical isochore % = °°. This 
implies that the expression for p(p, T) and/;(p, T) in 
the (jc, t) plane exist and are everywhere convergent 
for large R < x ^ oo. 

In general, the isothermal compressibility Kj is 
given by 






|Ap 



8-1 



I (Ap)' 



(8 + 5- 



1)-*-*- 

; dx 



fJL s (x) (16) 



with the thermodynamic stability [5] requiring that K T 
be everywhere non-negative. The dominant thermo- 
dynamic behavior near the critical isochore (x large) 
is given by the leading term in this expression (s= 1) 
and is like that found from simple scaling. As the 
critical point is approached along the coexistence 
curve within the single-phase region, the compress- 
ibility diverges as \t\~ y ' where y' = )3(8 — 1 ) , provided 

p\ ( — X ) 7^ 0. 

The thermodynamic behavior of the specific heat 
in the homogeneous phase is given by 



C v =-T(d 2 F/dT 2 )p, 
which may be written as 

p's(x)— PciMi'M+M'i'-iM} 



1 + f x 

__ j? |f|«+*-l>-2 
= 1 



PcTc£ 



M&+S-D- 



(17) 



with /JL(){x)=0 (neglecting the smooth contribution 
from F (p, 7 1 )). Since p\ {x) = p c p x (x) , the lowest 
possible contribution to the specific heat must come 
from an 5^2 term in this expression. Here in the 
single phase, just as in the two-phase region, the 
available experimental evidence suggests that the 
pressure term (d 2 pldT 2 )p c is much larger than the 
chemical potential term [12]. About the critical iso- 
chore above T c , the specific heat varies as |^|~ a , 
where a = 2 — (3(8 + s— 1), with s determined by the 
first nonvanishing term in the limit as %— >°°. (The 
anomaly in the specific heat along the critical isochore 
is also given by F r \t\~ a , where a=2 — j8(8+r) andF r 
denotes the limit of f r (x) ■ x -p(d+r-i) as jc— *°°.] Such 
structure admits the possibility of different values for 
the critical indices a and a'. If, for example, the 
leading s' = 2 were finite in the two-phase region and 
then vanish above T c at large x [see eqs 15, 17], s 



would be greater than 2 and thus the difference 

\a — a' | =/3|s — s'| could be greater than zero. 



3. Discussion 

The extended formulation of thermodynamic scaling 
presented here provides a quite general description 
for the critical region of a fluid. The formalism proceeds 
beyond the usual lowest order asymptotic behavior 
and accommodates both the chemical potential and 
pressure as equivalent physical variables. Within the 
description it is possible to note contributions from the 
various higher ordered terms {s > 1) and of the two 
quantities p, p to the various physical anomalies. We 
found, for example, that the difference between the 
potential and pressure variables is not only reflected 
in the shape of the coexistence curve (as determined 
by the value of Xi), but also enters into the leading 
behavior of the specific heat. 

We also expect extended scaling to be valid beyond 
the asymptotically small range of simple scaling, with 
contributions from the additional terms affecting the 
determination of the values of the critical indices away 
from the acutal critical point. For example, the experi- 
mentally determined shape of the coexistence curve 
/3* would be greater or less than the true /3, the modi- 
fication depending upon both the sign and magnitude 
of the additional terms of eq (1). 

In the rather special cases where 1//3 is an even 
integer, the formulation recovers the classical thermo- 
dynamic description of a van der Waals system. For 
j3 =1 /2, this is just the familiar parabolic coexistence 
curve of the mean-field result, with the usual recti- 
linear diameter p = p c + p->\t\. 

As mentioned before, any comparison of these re- 
sults to the physical features must remain somewhat 
qualitative. This is because the formulation is essen- 
tially analytic everywhere and must therefore suffer 
the same defects as all "classical" equations of state. 
However, we believe the utility of this analysis lies 
in the display of certain features which must eventually 
be contained in any complete thermodynamic descrip- 
tion of the critical region. 
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